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ABSTRACT 

This is the second paper of a series that reports on our investigation of the clustering properties of 
AGNs in the R OS A T All-Sky Survey (RASS) through cross-correlation functions (CCFs) with Sloan 
Digital Sky Survey (SDSS) galaxies. In this paper, we apply the Halo Occupation Distribution (HOD) 
model to the CCFs between the RASS Broad-line AGNs with SDSS Luminous Red Galaxies (LRGs) 
in the redshift range 0.16 < z < 0.36 that was calculated in paper I. In our HOD modeling approach, 
we use the known HOD of LRGs and constrain the HOD of the AGNs by a model fit to the CCF. 
For the first time, we are able to go beyond quoting merely a 'typical' AGN host halo mass, Mh, 
and model the full distribution function of AGN host dark matter halos. In addition, we are able to 
determine the large-scale bias and the mean Mh more accurately. We explore the behavior of three 
simple HOD models. Our first model (Model A) is a truncated power-law HOD model in which all 
AGNs arc satellites. With this model, we find an upper limit to the slope (a) of the AGN HOD that 
is far below unity. The other two models have a central component, which has a step function form, 
where the HOD is constant above a minimum mass, without (Model B) or with (Model C) an upper 
mass cutoff, in addition to the truncated power-law satellite component, similar to the HOD that is 
found for galaxies. In these two models we find that the upper limits on a are still below unity, with 
a < 0.95 and a < 0.84 for Model B and C respectively. Our analysis suggests that the satellite AGN 
occupation increases slower than, or may even decrease with, Mh in contrast to the satellite HODs of 
luminosity-threshold samples of galaxies, which, in contrast, grow approximately as (N s ) oc M^ with 
a ~ 1. These results are consistent with observations that the AGN fraction in groups and clusters 
decreases with richness. 

Subject headings: galaxies: active — X-rays: galaxies — cosmology: large-scale structure of Universe 



1. INTRODUCTION 

Investigating how active galactic nuclei (AGNs) are 
distributed in the Universe provides a clue as to the 
physical conditions in which accretion onto supermassive 
black holes (SMBH) takes place. After the formation of 
galaxies in the universe, at any given time in the his- 
tory of the Universe, only a small fraction of galaxies 
show AGN activity. It is known that almost all galax- 
ies with a spheroidal component contain a SMBH, and 
that the BH mass is closely related to the mass or the 
velocity dispersion of the spheroidal component (e.g., 
Ferrarese & Merritt 2000; Gebhardt et al. 2000; Mar- 
coni & Hunt 2003; Haring & Rix 2004; Shankar 2009 
for review). This means that most galaxies have likely 
had one or more brief AGN periods, during which the 
central SMBH grows in such a way that the growth of 
the SMBH and the formation of the spheroidal compo- 
nent are tightly related. The question of when and un- 
der what physical conditions the accretion takes place 
is important in understanding not only the origin and 
evolution of SMBHs but also the origin and evolution 
of galaxies. Many models postulate that AGN activity 
is merger-driven, especially those with high luminosities 
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(i.e., QSOs), (e.g., Wyithe & Loeb 2002; Hopkins et al. 
2006). On the other hand, for lower-luminosity AGNs, 
processes internal to the galaxy such as galaxy disk in- 
stabilities, may be important (Kauffmann et al. 2007; 
Hasinger 2008). Different mechanisms may be responsi- 
ble for triggering AGN activity at different redshifts or 
luminosities. 

Large-scale clustering properties of AGNs provide im- 
portant clues as to the physical process(es) that are re- 
sponsible for the SMBH accretion. The large-scale clus- 
tering amplitude reflects, through the bias parameter &, 
the typical mass of the dark matter halos (DMH) in 
which AGN reside (Sheth & Tormen 1999; Sheth et al. 
2001; Tinker et al. 2005). Thus exploring the clustering 
properties of AGNs at different redshifts, luminosities, 
and AGN types to probe the masses of the DMHs that 
host them, provides us with important clues in under- 
standing constraints on SMBH growth in a cosmological 
context. 

Clustering measurements of X-ray point sources have 
been made whenever new large-scale surveys have been 
completed, through angular correlation functions (e.g. 
Puccetti et al. 2006; Miyaji et al. 2007; Ueda et al. 
2008; Plionis et al. 2008; Ebrero et al. 2009b). When 
redshifts for a complete sample become available, three- 
dimensional (3D) correlation functions can then be com- 
puted (e.g., Mullis et al. 2004; Yang et al. 2006; Gilli et 
al. 2009; Cappelluti et al. 2010). Redshift information 
provides a huge leap in clustering measurements both in 
terms of the statistical accuracy as well as in removing 
systematic errors associated with model assumptions in 
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making Limber's de-projection (Limber 1954). However, 
even with redshift information, small number statistics 
has limited the accuracy of correlation function (CF) 
measurements, especially through the auto-correlation 
function (ACF) of the AGNs themselves. The situa- 
tion can be improved by measuring the cross-correlation 
function (CCF) of AGNs with a galaxy sample that has 
a much higher space density. In order to facilitate a 
CCF analysis, one needs an extensive galaxy redshift 
surveys with common sky and redshift coverage as the 
AGN redshift surveys. With recent large-scale survey 
projects, measurements of AGN clustering through CCF 
with galaxies are now emerging (e.g., Li et al. 2006; Coil 
et al. 2007, 2009; Hickox et al. 2009; Padmanabhan et al. 
2009; Mountrichas et al. 2009). 

While the ROSAT All Sky Survey (RASS) (Voges et 
al. 1999), most of which was conducted during the first 
half year of the ROSAT mission (mostly in 1990), pro- 
duced a catalog of « 120, 000 X-ray point sources, the 
availability of a comprehensive redshift survey of RASS- 
selected AGNs was limited until a RASS-Sloan Digital 
Sky Survey (SDSS) matched AGN catalog became avail- 
able (Anderson et al. 2003, 2007). In view of this, we 
initiated a series of studies of the clustering properties 
of low redshift AGNs in the RASS-SDSS sample through 
the CCF analysis with SDSS galaxies. In Krumpe et al. 
(2010, hereafter paper I), we reported our results on our 
CCF analysis between broad-line AGNs in RASS that 
have been identified with SDSS (Anderson et al. 2007) 
and the SDSS Luminous Red Galaxies (LRGs) (Eisen- 
stein et al. 2001) in the redshift range 0.16 < z < 0.36. 

Through these efforts to measure the AGN ACF/CCF, 
the masses of the DMHs where AGN activity occurs are 
being gradually uncovered. Most results measuring the 
3D correlation functions of X-ray selected AGNs indicate 
that the typical DMH mass in which these AGN reside 
is in the range 12 < log M h [/i _1 M sun ] < 13.5 at both 
low (z < 0.4) and high (z « 1) redshifts (e.g., paper 
I; Coil et al. 2009; Gilli et al. 2009; Cappelluti et al. 
2010). Optically-selected QSOs, typically representing a 
high redshift, high luminosity AGN population, are as- 
sociated with DMHs with a typical mass in the range 
12 < log M h [h~ 1 M sun ] < 13 (e.g. Porciani et al. 2004; 
Coil et al. 2007; Shen et al. 2007; Ross et al. 2009). It 
appears that the luminosity dependence of AGN clus- 
tering may be different at different redshifts. In paper 
I, we show that among broad-line AGN, high X-ray lu- 
minosity AGNs arc more strongly clustered than lower 
X-ray luminosity AGN at z« 0.3. At z ~ 1, however, 
low luminosity X-ray selected AGNs are more strongly 
clustered than optical QSOs (Coil et al. 2009). This ap- 
parent difference may be caused by a non-monotonic lu- 
minosity dependence and/or AGN type dependence of 
biasing rather than a dependence on redshift. What is 
needed is to explore AGN clustering across wider ranges 
in luminosity-redshift space to break this degeneracy. 

In interpreting the correlation function measurements, 
most previous studies use the large scale bias (b) of AGNs 
to infer the associated typical DMH mass using linear 
growth and linear biasing schemes. Strictly speaking, 
this is only valid for the correlation function measure- 
ments on sufficiently large scales (r > 1 — 2h~ 1 Mpc). 
Non-linear modeling through the Halo Occupation Dis- 
tribution (HOD) framework (e.g. Peacock & Smith 2000; 



Seljak 2000; Cooray & Sheth 2002) is imperative to accu- 
rately interpret and make full use of the correlation func- 
tion measurements. In this framework, the mean number 
of the sample objects in the DMH ((A)) is modeled as 
a function of the DMH mass (Mh). Then the two-point 
correlation function is modeled as the sum of the contri- 
butions of pairs from the same DMH (1-halo term) and 
those from different DMHs (2-halo term). This method 
has been used extensively to interpret galaxy correlation 
functions (e.g. Hamana et al. 2004; Tinker et al. 2005; 
Phleps et al. 2006; Zheng et al. 2007; Zehavi et al. 2010; 
Zheng et al. 2009, hereafter Z09) to constrain how various 
galaxy samples are distributed among DMHs as well as 
whether these galaxies occupy the centers of the DMHs 
or are satellite galaxies (Kravtsov et al. 2004; Zheng et 
al. 2005). 

Partially due to the low number density of AGNs, 
there have been few results in the literature interpret- 
ing AGN correlation functions using HOD modeling, 
where the small-scale clustering measurements are es- 
sential. Padmanabhan et al. (2009) discussed qualitative 
HOD constraints on their LRG-optical QSO CCF, where 
they argued that the acceptable models include those in 
which > 25% of their QSOs are satellites and those in 
which QSOs are a random subsampling of a luminosity- 
threshold sample of galaxies. Shen et al. (2010) also used 
the HOD modeling approach to binary pairs of QSOs at 
3 < z < 4 and conclude that they favor a model in which 
> 10% of the QSOs are satellites. 

In this paper, we apply the HOD model to the CCFs 
of the RASS broad-line AGNs (hereafter, simply referred 
to as AGNs) and SDSS LRGs obtained in paper I. The 
HOD modeling has allowed us to model the CCFs be- 
yond simple power-law fits and to investigate constraints 
on how these AGNs are distributed among DMHs as a 
function of halo mass. While almost all previous HOD 
modeling studies apply the method to ACFs, here we ap- 
ply the HOD modeling approach to the CCF. With this 
method, we use the previously estimated LRG HOD to 
model the AGN HOD by fitting to the measured AGN- 
LRG CCF. In our analysis, we take advantage of the LRG 
HOD of Z09, which is based on the well- measured LRG 
ACF from SDSS (Zehavi et al. 2005b). This approach 
was taken by Z09 in modeling LRG-L* galaxy CCFs, as- 
suming that the LRGs are central galaxies of the DMHs 
and that the L* galaxies are satellite galaxies upon the 
calculation of the 1-halo term. 

In this paper, we present a comprehensive explanation 
of the HOD analysis that can be applied to general cases, 
and then we apply the method to our AGN-LRG CCF. 

The scope of the paper is as follows. In sect. 2, we 
summarize the samples used in paper I and the basic 
methods of CCF calculations. In Sect. 3, we explain 
our basic modeling procedure of the LRG ACF, which is 
used as a template to apply the HOD modeling to the 
AGN-LRG CCF. In Sect. 4, we show our results on the 
AGN HODs. Sect. 5 discusses our results, including a 
comparison with the results from paper I and the as- 
trophysical implications of our constraints on the AGN 
HOD obtained here. Finally, Sect. 6 summarizes impor- 
tant consequences from our analysis and concludes our 
discussion. 

Throughout the paper, all distances are measured in 
comoving coordinates and given in units of /i _1 Mpc, 
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where h = _ffo/100kms _1 . We also use the symbol 
hro = -ffo/70kms _1 for X-ray luminosities, while h = 1 
is used to express optical absolute magnitudes, for consis- 
tency with referenced articles. We use a cosmology with 
h = 0.7, ft M = 0.3, tt A = 0.7, a 8 = 0.8, n haryon = 0.047 
and n s = 1, which are consistent with the most updated 
WMAP cosmology as of writing this paper (Spergel et 
al. 2003) 4 . The symbol log signifies a base-10 logarithm, 
while a natural logarithm is expressed by an In. 

2. THE RASS AGN-LRG CCF MEASUREMENTS 

In paper I, we calculate the CCFs between AGNs and 
the LRG, as well as the LRG ACF in the redshift range 
0.16 < z < 0.36. Basic properties of the sample used in 
that paper are repeated here in Table 1. The rationale 
and strategies of the sample selection are discussed in 
paper I in detail, including a redshift-_L x diagram of the 
AGN sample (Fig. 1 of paper I). Here we briefly sum- 
marize the sample selection and CF calculations below. 

We use the AGN sample from the RASS-SDSS 
matched AGN catalog by Anderson et al. (2003, 2007), 
which is based on the SDSS Data Release (DR) 5. In 
order to make a well-defined, uniform AGN sample, we 
only include broad lines AGNs. For the reference sam- 
ple, we extract LRGs from the SDSS Catalog Archive 
Server Jobs System 5 using the flag "galaxy _red" , which 
is based on the selection criteria defined in (Eiscnstein et 
al. 2001). We create a volume-limited spectroscopic LRG 
sample with -23.2 < M g z=0 - 3 < -21.2 and 0.16 < z < 
0.36. The AGN sample is also limited to this redshift 
range. In order to obtain the largest common geome- 
try between the AGNs and LRGs with public geometry 
and completeness files (Blanton et al. 2005), we limit our 
sample to the SDSS DR4+ geometry 6 . The selected ge- 
ometry covers an area of 5468 deg 2 . The AGN sample is 
further divided into high Lx and low Lx samples. The 
numbers of LRGs and AGNs used in our CCF measure- 
ments are listed in Table 1. 

In calculating the AGN-LRG CCFs and the LRG ACF, 
we use the classic Davis & Peebles (1983) estimator, in 
a two-dimensional (2D) grid in the (r p ,n) space, where 
r p is the projected distance and ir is the line-of-sight 
separation (both in comoving coordinates): 



DiD 2 (r p ,ir) 
DiR 2 (r p ,Tr) 



1, 



(1) 



gas in our galaxy, combined with the diversity of the 
X-ray spectra of AGN, makes it very difficult to model 
the spatial variation in the detection limit of the AGNs, 
which is essential in generating a random sample. 
We then calculate the projected correlation function: 



w P (r p ) 



£(r p ,7r)d7T, 



(2) 



where the upper bound of the integral (7r max ) is de- 
termined by the saturation of the integral. We take 
TTmax = 80 and 40/j- 1 Mpc for the LRG ACF and AGN- 
LRG CCF respectively. The reasoning for these choices 
is explained in paper I. 

In paper I, we further calculate the inferred AGN 
ACF, w p (AGN\AGN), from the AGN-LRG CCF 
w p (AGN\LRG) and LRG ACF w p (LRG\LRG), assum- 
ing a linear biasing scheme: 



w p (AGN\AGN) = 



[w p (AGN\LRG)Y 
w p (LRG\LRG) 



(3) 



in which the subscripts 1 and 2 identify the sample, 
DiD 2 is the number of pairs between real samples, and 
DiR2{r p ,Tr) is the number of pairs between real sample 
1 and random sample 2. For the AGN-LRG CCF, the 
sample 1 is the AGN sample and the sample 2 is the LRG 
sample, while for the LRG ACF, both sample 1 and sam- 
ple 2 are the same LRG sample. We use Eq. 1 rather 
than that proposed by Landy & Szalay (1993) because 
Eq. 1 requires a random sample for only sample 2 (in 
our case, LRGs). This is important due of the difficulty 
in generating a random sample for a RASS-based AGN 
catalog. ROSAT is sensitive to soft X-rays, which are 
subject to absorption by neutral hydrogen. The varia- 
tion in the column density across the sky due to neutral 

4 http: / /lambda. gsfc. nasa. gov / product / map / current / parameters . cfm 

5 http://casjobs.sdss.org/CasJobs/ 

6 http://sdss.physics.nyu.edu/lss/dr4plus 



We then fit w p (AGN\ AGN) with a power-law form and 
derive the AGN bias parameters and typical masses of 
the DMH in which the AGNs reside, based on the power- 
law fits. 

In this paper, instead of using the inferred 
AGN ACF as in Paper I, we directly model the 
CCF, w p (AGN\LRG), with an HOD ' analysis, as 
w p (AGN\LRG) has been directly derived from the ob- 
servations and does not depend on a linear biasing ap- 
proximation. The comparison of results from paper I and 
this paper are discussed below in Sect. 5.1. 

3. MODELING THE HALO OCCUPATION DISTRIBUTION 

3.1. Model Ingredients 

In performing the HOD modeling, we consider that 
galaxies and AGNs are associated with DMHs, the mass 
function of which (per comoving volume per dM^) is de- 
noted by 0(Mh) cLMh, where Mh is the dark matter halo 
mass. We use <j>(Mh) based on Sheth & Tormen (1999), 
which is in good agreement with later analyses by Sheth 
et al. (2001) and Jenkins et al. (2001). A DMH may con- 
tain one or more galaxies and/or AGNs that are included 
in our samples. We then model the two-point projected 
correlation function as the sum of two terms, 

w p (r p ) = w p ,ih(r p ) + w p ,2h(r p ), (4) 

where the terms are: 

• the 1-halo term (denoted by a subscript "Hi"), 
where both of the objects occupy the same DMH, 
and 

• the 2-halo term (denoted by a subscript "2h"), 
where each of the pair of objects occupies a dif- 
ferent DMH. 

Recent articles prefer to define £ih such that the total 
3D correlation function is expressed by £ = (l+£ih)+£2h 
(Tinker et al. 2005; Zheng et al. 2005; Blake et al. 2008), 
instead of £ = £ih + ^2h, as used in older articles. This is 
because 1 + £ represents a quantity that is proportional 
to the number of pairs, and allows one to express each 
of the 1- and 2-halo terms consistently. In this case, the 
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TABLE 1 

Properties of the LRG and RASS-AGN samples 



sample 
name 


z 

range 


Mg range 
[mag] 


number 


< "LRG > 

[h 3 Mpc- 3 ] 


< 2 > 


^ A/fZ — 0.3 \ 

[mag] 


LRG sample 


0.16 < z < 0.36 


-23.2 < M|= - 3 < -21.2 


45899 


9.6 x 10~ 5 


0.28 


-21.71 




z-range 


logLx range a > b 
range [/tf () 2 crg s _1 ] 


number 


< "AGN > 

[h 3 Mpc' 3 ] 


< z > 


log < L x > a ' c 
[h-m erg s" 1 ] 


All RASS-AGN sample 
High L x RASS-AGN sample 
Low L x RASS-AGN sample 


0.16 < z < 0.36 
0.16 < z < 0.36 
0.16 < 2 < 0.36 


43.7 < logLx < 44.7 
44.3 < logLx < 44.9 
43.6 < logLx < 44.3 


1552 
562 
990 


6.0 x 10~ 5 
1.2 x 10~ 6 
5.8 x 10~ 5 


0.25 
0.28 
0.24 


44.17 (44.16) 
44.58 (44.53) 
43.95 (44.16) 



a X-ray luminosity is measured in the rest-frame energy range 0.1-2.4 keV. 

b The symbol "<" is used in case of a "soft" luminosity boundary and indicates the 10th (for the lower 
bound) or the 90th (for the upper bound) percentile from the lowest luminosity object in the sample. The 
symbols < and < are used to indicate a hard luminosity boundary we impose. 

c The the logarithm of the mean Lx is followed by parenthesized median values. 



number of pairs is oc 1 + £ = [1 + £ih] + [1 + ^2h]- In 
this new convention, our w P! ih represents the projection 
of 1 + £ih rather than £ih, i.e., 



/oo 
[1+Cih (r p ,7r)]ti7T. 
-OO 



(5) 



Similarly, we express the power spectrum of the distri- 
bution of the objects in terms of the 1- and 2-halo term 
contributions: 



with 



P(k) = P lh (k) + P 2h (k), 



/k 
—Pih(k)J (kr p )dk 



(6) 



w Pl 2h(? , p)= / —P2h{k)Jo(kr p )dk, 



(7) 



where Jo(x) is the zeroth-order Bessel function of the 
first kind. 

In the case of modeling the ACF of a sample, the 2- 
halo term of the power spectrum can be approximated 
by (Cooray & Sheth 2002) 

b 2 P lin (k,z), 



(8) 



where b is the bias parameter of the sample, which we 
model as 



J b h (M h )(N){M h )(j ) (M h )dM h 
f(N)(M h )0(M h )dM h ■ 



(9) 



For the linear power spectrum, Pu n (k, z), we use the pri- 
mordial power spectrum with n s — 1 and a transfer func- 
tion calculated using the fitting formula of Eisenstein & 
Hu (1998) under our assumed cosmology (Sect. 1). For 
the mass-dependent bias parameter of DMHs, &h(Mh), 
we use Eq. (8) of Sheth ct al. (2001) with recalibrated 
parameter values by Tinker et al. (2005) (see their Ap- 
pendix A). 

In calculating the 1-halo term, the usual assumption 
is that the radial distribution of the involved objects 
that are not at the center of the halo - objects that are 
satellites - follows the mass profile of the DMH itself, 
i.e., there is no internal biasing. Following Zheng et al. 
(2009), whose results will be used later in this work, we 



use the Navarro et al. (1997) (NFW) profile for the DMH 
density distribution, which is still a popular choice for 
representing DMH profiles, although other parametriza- 
tions exist in the literature (e.g. Knollmann et al. 2008; 
Stadel et al. 2009). We express the Fourier transform of 
the NFW profile of the DMH with mass M^, normalized 
such that volume integral up to the virial radius is unity, 
by y(k,Mh). Again, in the case of the ACF, 



Plh(k): 



1 



<KM h )[(N c N s ){M h )y(k,M h ) + 



(2^) 3 n 2 

(N S (N S - l))(M h )\y(k,M h )\ 2 ]dM h , (10) 

where N c and N s are the numbers of the objects in 
the sample per DMH function of Mh for those 

that are at the center of the halos (central) and 
those that are off center (satellites) respectively, while 
n = J (N)(Mh)(j)(Mi 1 )dMi 1 is their overall space density 
(Hamana et al. 2004; Seljak 2000). 

The HOD model is calculated using a set of software 
we have developed, based partially on a code developed 
by J. Peacock for the use in Peacock & Smith (2000) and 
Phlcps et al. (2006). Our software includes a number of 
improvements over their code, most notably by including 
the application of the HOD modeling to a CCF. We note, 
however, that in our analysis we do not include some re- 
cent improvements in modeling the 2-halo term (Zheng 
et al. 2005; Tinker et al. 2005; Zheng et al. 2009). The 
neglected factors include the convolution of the DMH 
profile to the 2-halo term and the scale-dependent bias. 
In addition, if a pair of objects are closer than the sum 
of the virial radii of their parent halos, the pair should 
not be counted in the 2-halo term; this effect is neglected. 
These improvements have enabled those authors to make 
accurate modelings of, e.g., the ACF of SDSS galaxies, 
where the correlation functions can be measured to a few 
percent level at r p w lh" 1 Mpc. However, as AGNs have 
a much smaller number density, the errors on our CCF 
measurements are > 20% (paper I) at this r p , where 
these improvements are important. Therefore the effects 
of these recent improvements by other authors are es- 
timated to be within the la measurement errors of the 
CCF. 

3.2. The HOD modeling of the AGN-LRG CCF 
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In this subsection, we explain the application of the 
HOD modeling to the CCF between AGNs and LRGs (or 
more generally, between any two different populations). 
Using the HOD of the LRGs as a template, which has 
been accurately constrained by Z09 using the LRG ACF 
measured by Zehavi et al. (2005b), we constrain the HOD 
of the AGNs by fitting to the AGN-LRG CCF. A similar 
approach has been made by Z09 to model the CCF be- 
tween LRGs and L* galaxies, where they assumed that 
the LRGs are at the halo centers and the L* galax- 
ies are satellites in the 1-halo term calculation. Since 
LRGs occupy the centers of almost all massive DMHs 
(logMh > 14[/i _1 M Q ], Z09), this is a good approxima- 
tion for their LRG-L* galaxy CCF as well as for our 
AGN-LRG CCF. However, for completeness, we develop 
a formulation of the HOD modeling that includes more 
general cases. Hereafter, the quantities for AGNs are 
represented by a subscript "A" , LRGs by "G" (repre- 
senting galaxies), and CCF between the two by "AG". 
We denote the LRG HODs at the halo center and of 
satellites by (AT Gc )(Mh) and (JVg s )(Mi,), respectively. 
Now (N G ){M h ) = (JV G ,c>(M h ) + <jv G ,s>(M h ). Likewise, 
the HODs of the AGNs at the halo centers, of satel- 
lites, and the sum of the two are denoted by {N A c )(Mh), 
(A A , s )(M h ) and (N A )(M h ). 

In our approximate treatment, the two halo term can 
be expressed as follows. Let the linear bias parameters 
of the LRGs and AGNs be b G and 6a respectively: 
_ / b h (M h )(N G )(M h )^(M h )dM h 
G f(N G )(M h )d>(M h )dM h [ ' 

J b h (M h ){N A )(M h )<j>(M h )dM h 
bA J(N A )(M h )cj ) (M h )dM h ■ ( ' 

Then the two halo term of the power spectrum corre- 
sponding to the CCF (cross power spectrum) can be ex- 
pressed by 

PAoMk) w b A b G P lin (k)- (13) 
The 1-halo term is composed of three terms: 

/ 5 AG,ih(fc) = 7 ^ / 0(M h ) X 

[<A A ,cA G , s + A A , s A G , c )(M h ) y(k, M h ) + 

<A A , s A G , s )(M h ) |y(fc,M h )| 2 ]dM h , (14) 
provided that no object is in common between the AGN 
and LRG samples (which is the case in our CCF) . In this 
case, the relation 

<A A , x A G , y )(M h ) = (A A , x )(M h )(A G , y )(M h ) (15) 
holds, where the subscripts x and y represent any combi- 
nation of the subscripts s and c, except for the case both 
are c's. This relation is exact even if (N) < 1, where one 
has to use the exact Poisson distribution for the proba- 
bility distribution of having exactly N objects in a halo 
(sometimes called a sub-Poissonian case). The relation 
(A A , c A GiS + A A , s A GiC ) = (Aa, c A GiS ) + (A A , S A G , C ) also 
holds. 

Even in cases where there are objects in common be- 
tween the two samples, the above discussion does not 
change for the central-satellite pairs, because the two 
objects are surely different. However, the existence of 
the common objects makes assumptions behind the third 
term in the integrand of Eq. 14 invalid. A detailed dis- 
cussion of this case is beyond the scope of the present 
paper, and will be made in a future paper. 



3.2.1. The LRG HOD 

The HOD of the LRGs has been extensively studied 
by Z09 based on the ACF measurements of the LRGs by 
Zehavi ct al. (2005b) using SDSS Data Release (DR) 3. 
Their (Af GiS )(Mh) model has essentially five parameters, 
which arc the values at five Mh points, and they have 
spline-interpolation values between these points. They 
take an elaborate parametrization of (A^ G>c )(Mh), which 
is an integration of a luminosity-dependent smoothed 
step-function, where the value increases from (at low 
Mh) to 1 (at high Mh), with a transition following a 
luminosity-dependent error function. 

In paper I, we re-calculate the LRG ACF for the 
0.16 < z < 0.36, -23.2 < M g < -21.2 sample based on 
DR4+, using «46,000 LRGs, instead of «30,000 LRGs 
from DR3 used by Zehavi et al. (2005b). We obtain prac- 
tically the identical ACF to that calculated by Zehavi et 
al. (2005b) . Since we want to use exactly the same model 
and the same LRG sample between the LRG ACF and 
the LRG- AGN CCF, we derive the best-fit LRG HOD 
to our DR4+ sample using our software. However, in- 
stead of newly exploring the full parameter space, we 
take advantage of the analysis in Z09 to find the HODs 
that give the best fit to our data by tweaking their LRG 
HODs as follows. Fig. 1 (b) of Z09 shows the A X 2 = 4 
upper and lower bounds of the HODs separately for the 
central and satellite LRGs. For the satellites, we start 
with their A% 2 = 4 upper and lower bounds, denoted by 
N% s and N% s respectively. We then interpolate between 
these curves: 

(A G , s )(M h ) = (1 - f)(Nl s )(M h ) + /(A z u s )(M h ), (16) 

where / is a weight parameter in the linear interpolating 
procedure. 

For the central LRGs, we shift their central HOD 
(Az jC ) horizontally by d in logM^ 

(JV G , c >(M h ) = <Az, c )(Mhl0 d ) (17) 

We tweak A GiC horizontally because of the (AT G;C )(Mh) < 
1 constraint. In calculating (N s (N s -l))(M h ) of the LRG 
ACF based on the above HODs, we calculate the mean 
value using the exact Poisson statistics in the small num- 
ber case ((N s ) < 10.0), while we use (N S )({N S ) - 1) for 
larger values. As a whole, our model has two "tweak" 
parameters, / and d. 

We calculate a series of model ACFs at z = 0.28 (the 
average redshift of our sample, see Table 1) in a parame- 
ter space grid exhaustively over a rectangular area in the 
(/, d) space that is large enough to include all the region 
with A\ 2 < 5. The grid spacings are 0.03 and 0.003 for 
/ and d respectively. 

Then we fit the model to the data by minimizing x 2 , 
taking the correlation of errors and the total number den- 
sity of LRGs (riLR G ) into account: 

X 2 = £{K>(rw) - w^\r p ,)]M^ x 

ij 

Mr P j)-< dl (rw)]} + 
(n L RG-< G ) 2 /< RG , (18) 

where the quantities from the model are indicated by 
a superscript "mdl", M is the covariance matrix, and 
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<7„ LRG is the lcr error of the LRG number density. The co- 
variance matrix and cr nLRG are estimated using the jack- 
knife resampling method (see paper I) . During the mini- 
mization process, the models w™ dl (r Pi j) are calculated by 
interpolating from the four nearest grid points in the pa- 
rameter space computed above. The minimization is per- 
formed using the MINUIT package 7 distributed as a part 
of the CERN program library. Including the number den- 
sity term in y 2 gives important parameter constraints, 
not only on the overall normalization, but also the shape. 
This is because Ag,c, by definition, has the absolute max- 
imum value of unity (as there can not be more than one 
central galaxy in a halo), and for the LRGs, the value 
saturates at 1 for M h > 10 14 /i _1 M Q . With this con- 
straint, nLRG becomes sensitive to the combination of 
the relative normalizations between Nq^ c and Ng iB and 
the "cut-off" mass of Ng >c . For the LRG sample de- 
fined in paper I, n LRG = (9.56 ±0.13) x I0~ 5 h~ 3 Mpc~ 3 . 
In performing the fits, we use only the data points that 
are dominated by either the 1-halo or 2-halo terms, ex- 
cluding the transition region (0.46 < r p [/i _1 Mpc] < 2.8). 
In this transition region, our model ACF becomes inac- 
curate due to our approximations, especially in neglect- 
ing the halo-exclusion effect in the 2-halo term, as de- 
scribed above in Sect 3.1. We use the data points down 
to r p = 0.2 [/i _1 Mpc], which is slightly below the lower r p 
limit used in paper I. Our measurements to this scale are 
consistent with those measure by Masjedi et al. (2006). 

Our best fit model has 'tweak' parameter values of / = 
1.21(1.12; 1.29) and d = -0.087(-0.095; -0.083), where 
the 68% confidence ranges for two interesting parame- 
ters (Ax 2 < 2.3) are given in parentheses. Our HODs are 
slightly higher than those measured by Z09, probably be- 
cause we include a data point at r p = 0.2 [/i _1 Mpc] and 
we use a slightly different linear power spectrum. The 
best-fit model w p (r p ) is compared with the observation 
as well as the Z09 model in Fig. 1(a). The correspond- 
ing HODs ((iV G , s (M h )), (N G , c (M h )) and the total) are 
shown in Fig. 1(b), and confidence contours in the (d,f) 
space are shown in Fig. 1(c). 

3.2.2. The AGN HOD 

With the LRG HOD in hand, we come to our main pur- 
pose of obtaining constraints on the HOD of the AGNs 
in the RASS-SDSS survey. 

As a parametrized form of the AGN HOD, we first try 
a simple truncated power-law form, assuming that all the 
AGNs are satellites (Model A): 

(N AtS ) cx M£Q(M h - M cr ), (JV A ,c> - (19) 

where ©(a;) is the step function (= 1 at x > 0; = at 
x < 0), M CI is a critical DMH mass below which the HOD 
is zero, and a is a power-law slope of the HOD above 
M CI . While formally we are assuming that all AGNs are 
satellites, the HOD constraints obtained from this as- 
sumption can be applied fairly to the sum of the central 
and satellite AGNs in cases where the AGN activity in 
the central galaxy of high mass log M^h -1 Mq] > 13.5 
DMHs is suppressed. For simplicity of explanation, we 
visualize the HOD of the central LRG as a step func- 
tion with a step at log M\ ] \h~ 1 M ( 7 ) ] = 13.5, and refer 

7 http:/ /wwwasdoc. web. cern.ch/wwwasdoc/minuit/minmain. html 



10 4 r(a) LRG ACF 1 
: Our model (lh+2h) : 



1-halo term 




1 10 
r p [i]-'Mpc] 




Fig. 1. — (a) The ACF of the luminous red galaxy (LRG) sample 
and the best-fit HOD model (upper panel) with residuals (lower 
panel). The data points and error bars show our measurements 
of the LRG ACF. The circles show the data points used for the 
fit, while the crosses show those excluded. The dashed (green), 
long-dashed (green) and thick-dashed lines show the 1-halo term, 
2-halo term, and the total respectively. The dot-dashed (magenta) 
curve represents the best-fit model by Z09. Our simplified model 
fits well in the regimes where the 1-halo term and 2-halo terms 
dominate, (b) Our best fit HODs for the central LRGs (red dashed 
line), satellite LRGs (blue long-dashed line), and the total, (c) 
Confidence contours in the d-f space (see text for details). The 
confidence levels are A% 2 = 1.0, 2.3 and 4.6. The cross shows the 
best-fit values. 

the DMH mass roughly above (below) this cut as higher 
(lower) mass. There is no contribution to the 1-halo term 
at lower mass DMHs, since the 1-halo term requires both 
AGNs and LRGs to reside at the same DMH mass. Since 
there is in practice no distinction between central and 
satellite AGNs in the 2-halo term (Eq. 13), whether the 
AGNs in low mass DMHs arc satellites or central makes 
no difference to the CCF, and indeed this model can 
fairly accommodate a case where a low-mass halo repre- 
sents one galaxy (halo) at the center. At higher masses, 
among the three terms in Eq. 14 the satellite AGN- 
central LRG pairs dominate the 1-halo term, because, as 
we show later, the present HOD analysis is sensitive at 
log M\ 1 [h~ 1 Mo\ < 14.0, where there are in practice no 
satellite LRGs (Fig. 1), and therefore the central AGN- 
satellite LRG term as well as the satellite AGN-satellite 
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LRG term are negligible. Thus Eq. 19 represents the 
sum of the central and satellite AGN HOD if there is no 
central AGN at higher masses. 

The underlying assumption in Model A is partially mo- 
tivated by observations of AGN "downsizing" (Ueda et 
al. 2003; Hasinger 2008; Ebrero et al. 2009a; Yencho et 
al. 2009), where the number density of high luminosity 
AGNs peaks earlier in the history of the universe and 
drops rapidly towards low redshift, while lower luminos- 
ity AGN activity peaks later. (We note, however, that 
a recent work by Aird et al. (2010) reports that this 
trend might be weaker than those reported previously.) 
One possible implication of this trend is that black holes 
at the centers of massive galaxies occupying the centers 
of (higher mass) DMHs stopped accreting long before 
z ~ 0.3, while accretion occurs more frequently in lower 
mass galaxies in satellites or at the centers of lower mass 
halos in the redshift range of our sample. Thus model A 
can be a demonstrative case of this scenario. Addition- 
ally, central galaxies of massive halos are mostly early- 
type luminous galaxies, where AGN activity is reported 
to be suppressed (Schawinski et al. 2010). 

In Sect. 4.3 we consider in detail other models (Models 
B and C) which uses well-studied galaxy HODs (central 
plus satellites) as templates and treat the separate cen- 
tral and satellite HODs explicitly. 

Using a fixed set of (N G , c )(M h ) and (N GjS )(M h ) de- 
rived in the previous subsection and the formulations 
developed earlier in this section, we calculate the ex- 
pected w Pj ag(^p) m the two parameter (M cr ,a) model 
of (NA){Mh)- Due to much smaller errors of the LRG 
ACF compared with those in the AGN-LRG CCF, we 
fixed the tweak parameter / and d at the best-fit values 
during the fits to the CCF. Shifting the tweak parame- 
ters to any point on the Ax 2 = 2.3 contour in Fig. 1(c) 
causes a shift of only A% 2 < 0.05 to the same AGN HOD 
model. This justifies the use of fixed LRG HODs during 
the fit to the CCF. 

We calculate a series of model CCFs in a parameter 
grid at z = 0.28 with the spacings of 0.1 and 0.05 for 
log M cr and a, respectively, over the rectangular region 
that we are interested in (see confidence contours in the 
next section) . While the mean redshift of the AGN sam- 
ples range from z — 0.24 to z = 0.28, the 1-halo and 
2-halo terms vary by only w 2% and w 0.1% respectively 
between these redshifts, justifying our single redshift cal- 
culations. We search for the best-fit model by minimizing 
the correlated \ 2 : 

X 2 = E{K»AG(rw) -^(rp.olMy 1 x 

ij 

[^(^•)-<ag(^)]}> (20) 

where the errors and covariance matrix are estimated us- 
ing the jackknife resampling method (paper I). Unlike in 
the case of LRGs, here we do not include the number den- 
sity term in determining \ 2 , because the CCF depends 
only on a and M cr in Eq. 19 but not the normalization. 
Thus the CCF constraints on a and M cr and the den- 
sity constraints on the normalization can be separated. 
After the constraints of these two parameters are made, 
the normalization can be determined by matching to the 
number density of the AGNs. We use only data points 
with r p > 0.3/i _1 Mpc, as smaller scale bins contain only 



a small number of AGN-LRG pairs (< 15 pairs for the 
total and much less in the high and low Lx RASS- AGN 
samples), which would not warrant the applicability of 
the x 2 statistic. Unlike in the case of the LRG ACF, 
we use all the data points in 0.3 < r p [h~ 1 Mpc] < 40. 
Due to the smaller mass for the DMH occupied by the 
AGN compared to satellite LRGs, the transition between 
1-halo and 2-halo dominated regimes is very narrow for 
the AGN-LRG CCFs, as will be seen in the next sec- 
tion. Even within this small transition region, the effects 
of the major source of inaccuracy in our approximation, 
i.e., halo-exclusion in the 2-halo term, is estimated be 
a few times smaller than our CCF measurement errors, 
based on the comparison of our LRG 2-halo term and 
that calculated by Z09. 

4. RESULTS FROM THE HOD MODELING 

4.1. Fit Results 

Using the methods described above, we obtain con- 
straints on the AGN HOD parameters log M cr and a in 
Eq. 19. For each point in this two parameter space, 
we can also calculate the mean dark matter halo mass 
occupied by the AGN sample, 

_ J M h (N A )(M h )<P(M h )dM h 
[ h> f{N A )(M h )<f>(M h )dM h ' 1 > 

as well as the effective bias parameter of AGNs 6a (Eq. 
12). 

Figure 2(a) compares our AGN-LRG CCF for the total 
RASS- AGN sample and our best-fit model. 

Figure 2(b) shows the confidence contours in the 
log M cr - a space, overlaid on underlying thin (green) 
contours showing the mean halo mass (Mh). A line of 
constant bias, 6a, roughly follows a constant (Mh) con- 
tour, with small deviations caused by the non-linearity of 
the &h(Mh) function. The degree of the deviation is such 
that log(Mh) for a given b decreases typically by « 0.4 
when a is decreased from w 1 to w — 3. Figures 3 and 4 
show our HOD fits and confidence contours, respectively, 
for the high Lx and low Lx RASS-AGN samples, with 
an emphasis on the comparison between the two. 

Figures 2(b) and 4 show that the parameter pair 
(log M cr , a) is tightly constrained roughly along the con- 
stant (Mh) line. This primary constraint comes from the 
amplitude of the 2-halo term, which is proportional to 
the bias parameter of the AGN sample (for the fixed bias 
of the LRG sample). In the HOD analysis, however, the 
inclusion of the 1-halo term adds additional constraints 
in the two-parameter space. 

In any of the total, high L x and low L x RASS-AGN 
samples, models where the number occupation of AGNs 
is proportional to Mh above M cr (a = 1) are excluded, 
while a constant number of AGNs per halo above M cr is 
preferred. Figs 2(b) and 4 put constraints on the values 
of M cr for each sample. Using the Ax 2 = 2.3 contour 
(68% confidence level for two parameters), the minimum 
mass is constrained to be log M cr [/i _1 M Q ] > 11.9, 11.9, 
& 10.7 for the total, high L x and low L x RASS-AGN 
samples respectively. For the total RASS-AGN sample, 
the contour is constrained to have a > —3.0, while for the 
high Lx and low Lx RASS-AGN samples, the contour 
allows the smallest a that we have explored (a = —3.25). 
This essentially gives constraints on the width of the 
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Fig. 2.— (a) The AGN-LRG CCF (for the total RASS-AGN sample) and our best-fit HOD model for Model A. The symbols and line 
styles are the same as in Fig. 1(a). (b) The confidence contours of the Model A parameters (log M cr ,a). The confidence contours (thick 
contours) are at the Ax 2 = 1.0, 2.3 and 4.6 levels. The thin-contours (green) show the mean halo mass log (M|,[/i _1 Mq]) derived from the 
model parameters. Every other contour level is labeled. The contours show that the parameter constraints are almost along a constant 
log {M^) line. 
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Fig. 3.— The AGN-LRG CCFs and the best-fit HOD models 
(Model A) for the (a) high L x and (b) low L x RASS-AGN samples. 
The meanings of line styles and symbols are the same as those in 
Fig. 2(a). 



HOD, and at least for the total RASS-AGN sample, a 
delta-function type HOD could be marginally excluded. 

The best fit HOD parameters, the mean DMH mass 
occupied by the AGNs, and the linear bias from the 
fits are summarized in Table 2. For the fitting parame- 
ters log M cr and a, we show the full range correspond- 
ing to 68% confidence range in the 2D parameter space 
(A% 2 = 2.3). As shown in Figs. 2 and 4, these two 
parameters are highly correlated and the confidence con- 
tours are skewed, thus one should be cautious in inter- 
preting these ranges. On the other hand, if we project 
the probability distribution on the variable (M n ) or the 
variable b (bias) , both of which are unique functions of 
our two fitting parameters, the projected probability dis- 



FlG. 4. — The confidence contours of the fits to the CCFs to 
Model A for the high Lx (thick blue lines) and the low Lx (thick 
red lines) RASS-AGN samples. The contours are at A\ 2 = 1.0, 2.3 
and 4.6 levels. The confidence contours are shaded for the ease 
of the visibility of the difference between the two samples. The 
log (M^) contours are shown in the same way as in Fig. 2(b). 

tribution becomes roughly Gaussian. Thus we take the 
A\ 2 < 1 range (68% error for one parameter) to esti- 
mate the la errors of these derived quantities, (Af n ) and 
6 A , in Table 2. 

We note that because of the exponential drop of the 
DMH mass function the HOD at very high M n con- 
tributes little to the CCF. In order to check the maxi- 
mum Mh that our analysis can reasonably constrain, we 
recalculate the CCF with the truncated power-law model 
with an upper mass cutoff: 

(JV A ) oc M h Q 6(M h - M cr )6(M cut - M h ) (22) 

We recalculate x 2 with varying M cut while fixing (M cr ,a) 
at the best-fit value above. The change of \ 2 from the 
M cut = oo is larger than unity at log M\ 1 [h^ 1 AI Q ] < 
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TABLE 2 

The HOD parameters for RASS AGNs at z ks 0.28 (Model A) 



RASS-AGN 
Sample 


logM cr a 
[fc^M©] 


a a 


(N A )(M CT ) h 


log <A/ h ) c 
r'M ] 


b A c 


x 2c 


Total 
High L x 
Low Lx 


12.6 (11.9; 13.0) 
12.9 (11.9; 13.2) 
12.5 (10.7; 13.0) 


-0.3 (0.4;-3.0) 
-0.5 (0.7;-3.25 d ) 
-0.3 (0.7;-3.25 d ) 


8xl0~ 2 
3.3xl0~ 3 
5.5xl0~ 2 


13.09±0.08 

12 Q7+0-17 


1.30±0.09 

I 44+0-04 
1.22±0.15 


(5.9 
10.9 
12.5 



a The 68% confidence range in the 2D parameter space (A^ 2 < 2.3). 
b The normalization at the nominal case. 
c The 68% error (A% 2 < 1) for one parameter. 
d The bound truncated at our parameter grid limit. 
6 The x 2 value calculated over 11 data points with a covariance matrix 

o 

-l 

-2 

-3 



14.0. Thus our HOD analysis is sensitive up to Mh 
10 14 /i _1 M Q , i.e., the mass of a poor cluster. 



4.2. HODs and Space Densities 

In order to illustrate the range of acceptable HODs, 
we plot a number of representative AGN HOD models 
accepted by our fits ((Na)(M\ 1 )) as function of DMH 
mass in Fig. 5(a). In Fig. 5(b), we also show the same 
sets of models in units of the spatial density per comoving 
volume per log of the DMH mass: 



<JV A )(M h ) 



d\og Mh 



(jV A )(M h )-ln(lO)M h 0(M h ). (23) 



These models have been normalized to the observed num- 
ber densities (see Table 1). In each panel, three curves 
are plotted, representing the best-fit case and two ex- 
treme cases. These are (1) for the best-fit values (solid 
lines), (2) the point on the A\ 2 = 2.3 contour that has 
the smallest M cr value (i.e., the upper left tip of the two- 
parameter 68% confidence area in Fig. 2(b) or 4) (dotted 
lines), and (3) the smallest a point on the A\ 2 = 2.3 (the 
bottom of the confidence region and narrowest possible 
HOD distribution, dashed lines). The best- fit model with 
a = —3.25 (the smallest a in our search grid) is used 
for the high Lx and low Lx RASS-AGN samples, for 
which the Ax 2 = 2.3 contour continues below a = —3.25. 
These three models have been plotted to illustrate the al- 
most full range of possible HODs statistically accepted 
by our analysis. In order to illustrate the DHM mass 
range that can be occupied by both LRGs and AGNs, 
and therefore contributes to the 1-halo term of the CCF, 
we also show the HODs of the LRGs for reference. Of 
course, these plots are for our rather restrictive trun- 
cated power-law model and by no means are intended to 
illustrate an exhaustive set of possible HODs. However, 
this illustrates that these three models have roughly the 
same average Mh, with varied widths of the HOD. An 
important constraint is that the AGN HOD can be no 
wider than the dotted line in 5(a)(b), which represents 
the point corresponding to the lowest extreme of M cr 
and the roughly highest extreme of a, among acceptable 
models with A% 2 < 2.3. The meaning of this constraint 
is discussed in Sect. 5.2. 

4.3. Effects of Central AGNs 

In this section we explore models in which central 
AGNs are explicitly included, using galaxy HODs as a 
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Fig. 5.— (a) The HODs {(N A ){M h )) of Model A for three 
representative points in the parameter space that are acceptable 
(Ax 2 < 2.3) from our fits for each of the total (top panel), high 
Lx (upper middle panel) and low Lx (lower middle panel) AGN 
samples. The solid, dotted, and dashed lines represent the (1) best- 
fit, (2) lowest M cr extreme, and (3) lowest a extreme respectively 
(see text). For reference, the bottom panel shows the LRG HODs 
for the central galaxies (short-dashed line), satellites (long-dashed 
line), and their sum (solid line), (b) The same set of models are 
shown in terms of the spatial number density per log M^. The 
meaning of the line styles is the same as those in panel (a). 

template. Galaxy HODs have been extensively investi- 
gated using, e.g., SDSS data, with a high statistical ac- 
curacy. A simple model for luminosity-threshold galaxy 
samples (i.e. more luminous than a given limit), intro- 
duced by (Zehavi et al. 2005a) is a three parameter model 
including a step function for the central HOD and a 
truncated power-law satellite HOD. We express the AGN 
HODs here with the same form, except that the global 
normalization is also a free parameter (Model B): 



(Aa,cH/a©(A4 

(^A,s) = / A 0(M h 



M min ), 

M min ) (Mh /Mi ) Qs , (24) 

where /a is the AGN fraction (duty cycle) among central 
galaxies at Mh > M m - m . We use the symbol a s to em- 
phasize that it represents the HOD slope for satellites. 
Mi is the DMH mass at which the number of central 
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Fig. 6. — (a) The confidence contours of the Model B param- 
eters in the (log M m j n ,a s ) space for Mi/M m i n = 23. (b) The 
confidence contours of the same set of parameters for Model C, for 
logM cc [/i _1 M0] = 13.0 and Mi/M min = 23. The contours arc at 
Ax 2 = 1.0, 2.3 and 4.6 levels. 

AGNs is equal to that of satellite AGNs. For luminosity- 
threshold galaxy samples, it is usually assumed that the 
central galaxy HOD saturates at unity at the high M n 
end (i.e., /a = 1 in Eq. 24), because central galaxies 
in the highest mass DMHs are expected to be luminous 
enough to be included in the sample (e.g. Sect. 3.2.1). 
Since only a small fraction of galaxies contain an AGN, 
this has to be multiplied by a factor /a(<< 1), the value 
of which does not affect the CCF as discussed above. The 
value of /a can be determined by normalizing the HOD 

to 77AGN- 

It is important to note the underlying physical assump- 
tion of Eq. 24, which is that the AGN duty cycle among 
central galaxies does not depend on the DMH mass. We 
explore the consequences of this model as an illustrative 
contrast to Model A, which assumes a suppression of the 
AGN duty cycle in central galaxies of high mass DMHs. 

We investigate the behavior of x 2 (Eq. 20) in the space 
of M m i n , Mi, and a s . Due to the low signal-to-noise ra- 
tio of the CCFs for our luminosity-divided subsamples, 
we limit our discussion to the total sample. Even with 
the total sample, our statistics do not allow us to con- 
strain all three parameters; therefore we investigate the 
constraints on two parameters by fixing the remaining 
parameter. The HOD analysis of luminosity-threshold 
samples with -21.5 < M r max < -21.0, studied by Zc- 
havi et al. (2005a) found M 1 /M min w 23 and a s « 1.2. 
We consider the typical HOD of SDSS galaxy samples in 
this luminosity threshold range as a template because 
their M min range (12.0 < log M min [/i- 1 M ] < 12.7) 
roughly coincides with the minimum DMH range for our 
AGNs, as we see from the range of M cr in the confi- 
dence contour in Fig. 2(b), which approximately coin- 
cide with the range of M m i n as we see below. We there- 
fore (i) fix Mi/M m - m = 23 and obtain constraints in the 
(log M m i n , a s )-space, and (ii) fix a = 1.2 and obtain con- 
straints in the (log M m ; n , Mi/M m ; n )-space. 

The resulting confidence contours for (i) are shown in 
Fig. 6(a). The best-fit parameters and x 2 values are 
summarized in Table 3 for both (i) and (ii), where we 
also show the a s = 1.0 case for (ii), in which the number 
of satellites is proportional to the DMH mass. 

For the Mi/M m - m = 23 case, the best-fit slope for the 



satellite is below unity (a s = 0.55) and x 2 = 6.9 is 
identical to that of Model A. The best-fit model with 
a s = 1.0 is marginally rejected (Ax 2 = 3.6) and the 
model with a s — 1.2 is significantly rejected (Ax 2 = 18). 
Models that are consistent with a s — 1.2 can be found 
at Mi/M m ; n > 30, but they are not consistent with 
the HOD of luminosity threshold galaxy samples in the 
-21.5 < M r max < -21.0 range. 

We also explore a model (Model C) where only lower 
mass DMHs contain a central AGN, where the physical 
picture is explained at the introduction of Model A (Sect. 
3.2.2). In this model, the central HOD in Model B is 
replaced by 

(N A ,c) = / A 0(M h - M min )0(M cc - M h ), (25) 

while the form of (Na, c ) stays the same. As a repre- 
sentative case, we choose log M cc [/i _1 M ] = 13.0, be- 
low which the DMH center is more frequently occupied 
by blue galaxies than red galaxies in the HOD analy- 
sis of color-separated samples by Zehavi et al. (2005a), 
while we keep Mi/M m i n = 23 (but see also Zehavi et al. 
(2010), which use a different parametrization on color- 
separated samples). The confidence contours are shown 
in Fig. 6(b) and the best- fit parameters and x 2 values 
are shown in Table 3. In this scenario, the best-fit slope 
for the satellite is a s — 0, while a s = 1.0 is rejected at 
the Ax 2 = 4.8 level. 

5. DISCUSSION 

5.1. Comparison with Results from Paper I 

Table 4 compares our results with those from paper I. 
The meanings of the table columns are: (1) sample, (2) 
the bias parameter from the HOD analysis (Eq. 12), (3) 
the bias parameter from paper I, i.e. for the power-law 
fit in 0.3 < r p [/i _1 Mpc] < 15, (4) the bias parameter 
for the power-law fit in 1.5 < r p [/j _1 Mpc] < 15 calcu- 
lated in the same way as paper I (see below), (5) the 
"mean" DMH mass from Eq. 21, (6) the "typical" DMH 
mass calculated from &a from the HOD model, defined 
by bh(M t yp) = 6a calculated with the bh(Mh) by Tinker 
et al. (2005) (see Sect. 3.1). (7) the "typical" DMH mass 
from paper I (i.e. from power-law fits) calculated with 
the & h (M h ) by Sheth et al. (2001), and (8) the "typical" 
DMH mass from the paper I bias parameter re-calculated 
using the b h (M h ) relation by Tinker et al. (2005). While 
the results from the the HOD analysis presented here 
and the simple analysis of paper I agree well for the high 
Lx RASS-AGN sample, the HOD analysis gives larger 
bias parameters than those derived in paper I for the to- 
tal and low L x RASS-AGN samples at the level of 1.5 
a. Since these bias parameters have been derived from 
the same dataset with different methods, the discrepan- 
cies are systematic rather than statistical. One of the 
major findings of paper I is the detection of the X-ray 
luminosity dependence of the correlation function. The 
difference in correlation lengths (ro)> between the high 
Lx and low Lx RASS-AGN samples measured in paper 
I is w 2.5(7, based on power-law fits with a fixed slope of 
7 = 1.9. The differences of the AGN bias parameters &a 
and associated typical DMH mass -M^typ derived from 
power-law fits are « 1.8<r. In the HOD analysis, the sig- 
nificance of the difference in &a decreases to 1.1a, because 
it is constrained mainly by measurements at large scales 
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TABLE 3 

The HOD parameters for the RASS AGNs (Models B & C) 



Model 


RASS-AGN 
Sample 


logM min a 

[r'M ] 


M cc 


logMl/M min a 


a s a 


/A c 


x 2d 


B 


Total 


12.5 (12.2;13.0) 




1.36 (fixed) 


0.55 (0.95;-1.50 b ) 


2.9xl0~ 2 


6.9 


B 


Total 


12.4 (11.9;12.6) 




1.65 (1.51;1.83) 


1.20 (fixed) 


3.3xl0~ 2 


8.4 


B 


Total 


12.3 (12.1;12.7) 




1.55 (1.36;1.81) 


1.00 (fixed) 


2.6X10" 2 


7.7 


C 


Total 


12.6 (12.2;13.0 b ) 


13.0 (fixed) 


1.36 (fixed) 


0.0 (0.84;-1.50 b ) 


3.7xl0~ 2 


6.8 



a The 68% confidence range in the 2D parameter space (A% 2 < 2.3). 
b The bound truncated at our parameter grid limit. 

c The AGN fraction of the central galaxy evaluated at the best- fit parameters. 
d The x 2 value calculated over 11 data points with a covariance matrix. 



TABLE 4 
Comparion of Results from Paper I 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


RASS-AGN 








log(M h ) 


log(M h , typ ) a 


log(M h>typ ) a 


log(M h _ typ ) a 


Sample 


HOD 


PL 


PL,> 1.5 b 


HOD,tin05 c 


HOD,tin05 c 


PL,smt01 c 


PL,tin05 c 












[h- 


-M j 





Total 1.32±0.08 l.ll+g;^ 1.52±0.21 13.09±0.08 13.00±0.12 12.58±g;|| 12.701°,; 
High L x l-44±°-°* l-44t°|2 1.56 ±0.23 13.261°;^ 13.19l°;°j 13.10l° :2 j 13.271°;; 
Low L x 1.23±0.15 0.88l°;^ ^-o.l 12 - 97 -o.ll 12 - 87 -o.ll 11.83l°^ 55 12.031°;; 



a The symbol M typ represent the DMH mass satisfying bh(M typ ) = 6a- 

b The range of r p [/i _1 Mpc] used in the power-law limited to 1.5 < r p [ft _1 Mpc] < 15. 

c Conversion from bias to M h : smtOl (Sheth et al. 2001) or tin05 (Tinker et al. 2005). 



where the 2-halo term dominates (see below). However, 
in terms of M^typ derived from the HOD analysis, the 
difference is « 1.8a because of the additional constraints 
from the 1-halo term. 

There are a number of differences in the processes 
used to derive the bias parameter between paper I and 
this work. In paper I, the bias parameters are based 
on the power-law fits to Eq. 3 over scales of 0.3 < 
r p [/i _1 Mpc] < 15, assuming a linear bias over the en- 
tire scale range. Then the power-law models are con- 
verted into the rms fluctuations over 8 ft -1 Mpc spheres 
(c8,agn), and the bias parameter is calculated using 
o'8,agn/[o'o£ ) (z)], where D(z) is the linear growth fac- 
tor. The main source of discrepancy is that, while paper 
I also uses data points down to a scale of r p = 0.3ft _1 
Mpc to obtain 6, in the HOD analysis presented here the 
constraints on b mainly come from the 2-halo term dom- 
inated region (r p > 1.5/i _1 ), while smaller scale data 
constrain other variables. Fig. 3 shows that the dif- 
ference between the high Lx and low Lx RASS-AGN 
samples is more prominent in the 1-halo term dominated 
regime, at r p < 1ft, -1 , than in the 2-halo term dominated 
regime. This means that, for the AGN-LRG CCF, the 
1-halo term is more sensitive to differences in the AGN 
HOD. The reason for this is that while bh(Mh) increases 
only weakly with increasing with Mh, the 1-halo term is 
directly proportional to the number of AGNs that are in 
the same DMHs as those occupied by LRGs, i.e., very 
massive ones (logMh ^ 13.5[ft _1 M ]), i.e. the 1-halo 
term is more sensitive to the existence of AGNs in high 
mass halos. 



In the literature, relative or absolute bias parameters 
are often calculated based on results from power-law fits 
(e.g., paper I; Mullis et al. 2004; Yang et al. 2006; Miyaji 
et al. 2007; Coil et al. 2009). The power-law models of 
£(r) are usually converted to the rms fluctuation over 
8ft." 1 Mpc spheres or are averaged up to the distance of 
20ft,- 1 Mpc (£(< 20ft~ 1 Mpc)). While some authors use 
only large scales (r p > 1 — 2ft~ 1 Mpc) to ensure that only 
the linear regime is used, others include smaller scales. 
Usually the rationale for including smaller scales is to 
have better statistics, combined with the empirical ob- 
servation that the correlation functions are well fit with 
a power- law model over the scale of 0.1 — 0.3 < r p < 
10 - 20ft~ 1 Mpc (0.3 < r p < 15ft," 1 Mpc in the case of 
paper I). As a check, we have re-calculated the bias pa- 
rameter for power-law fits to the w p (AGN\ AGN) used in 
Paper I on scales of 1.5 < r p < 15ft -1 Mpc, to limit our- 
selves to the linear regime. These are shown in column 
(4) of Table 4. The bias parameters obtained in this man- 
ner carry large statistical errors, but the values are con- 
sistent with those derived from the HOD analysis within 
lcr in all three samples, while they are more than ss la 
above the bias parameters derived in paper I for the total 
and low L x RASS-AGN samples. The AGN bias param- 
eters in our total and high L x samples from the HOD 
analysis, as well as the power-law fit 1.5 < r p < 15ft _1 
Mpc, are consistent within combined la errors with those 
measured by Padmanabhan et al. (2009) for their QSO 
samples ALL0 and LSTAR0, which covers a similar red- 
shift range as ours (0.25 < z < 0.35). However, their 
bias parameters for their higher redshift QSO samples, 



12 



Miyaji et al. 



ALL1 & LSTAR1 (0.33 < z < 0.50) as well as ALL2 & 
LSTAR2 (0.45 < z < 0.60) arc lower. 

Using the power-law fit results including the non-linear 
regime is useful to detect the different clustering proper- 
ties between, e.g., different source populations, luminosi- 
ties, or types. In our case, this gives the most statistically 
significant difference between the clustering properties of 
the high Lx and low Lx RASS-AGN samples. However, 
one has to be careful in interpreting the results, e.g., in 
terms of the DMH mass based on the bias parameter. In 
principle, the HOD analysis overcomes this limitation, 
because we are able to derive the pure linear bias pa- 
rameter, while at the same time utilizing the data in the 
non-linear regime. 

5.2. Implications for the Environments of AGN 
Accretion 

With the HOD analysis, we can also obtain constraints 
on at least one additional variable which is indepen- 
dent from the "mean" DMH mass, thanks to the con- 
straints from the 1-halo term. The importance of the 
1-halo term constraints has also been noted by Padman- 
abhan et al. (2009) using a similar CCF analysis between 
photometric-redshift LRG samples and optically-selected 
QSOs in SDSS at z ~ 0.5. They qualitatively con- 
cluded that at least some QSOs must be satellites and 
also that models of the form of our Model B (Eq. 24) 
with Mi/M m ; n = 20 and a s — 0.9 are also acceptable, 
while lower a values are also possible. Our analysis here 
gives qualitatively similar results. 

The acceptable parameter space for Model A is always 
in the a < 1 regime for all of our three AGN samples: the 
total, high Lx, and low Lx RASS-AGN samples. This 
means that the average number of AGNs in a DMH docs 
not increase proportionally with mass. Figs. 2(b) and 4 
show that models with a < are preferred. It is interest- 
ing to compare this result with HOD analyses of galaxies, 
that find a s « 1 — 1.2 for a wide range of absolute mag- 
nitudes and redshifts at least up to z ~ 1.2 (Zehavi et al. 
2005a; Zheng et al. 2007; Zehavi et al. 2010). In princi- 
ple, neither the results nor the assumptions of Model A 
say anything about the central fraction among AGNs at 
the lower mass DMHs (see Sect. 3.2.2), and thus they 
do not directly constrain the slope of the satellite AGN 
HOD a s . However, if we take the extreme case where 
all the AGNs are satellites, a s < 0.4, this would imply 
that AGN fraction among satellite galaxies strongly de- 
crease with Mh- In Model B, where the central galaxy 
component is included explicitly under the assumption 
that the AGN duty cycle in central galaxies is constant 
against varying M h above M min , the best-fit slope is still 
less than unity (a s = 0.55), and models with a s w 1 
are only marginally rejected for Mi /M m ; n = 23 (a repre- 
sentative value of luminosity-threshold galaxy samples). 
In case of Model C, where only log M h [h~ 1 M Q ] < 13.0 
DMHs can contain a central AGN, the best-fit solution 
is a s = 0, while a s « 1 is rejected at a ~ 2a level 
for M\/M m i n = 23. On the whole, our results tend to 
favor models in which the AGN fraction among satel- 
lite galaxies decreases with increasing DMH mass (or 
the richness of groups and clusters), in the mass range 
12 - 13 < log Mh^Mg] < 14.0. 

Since our AGN sample is selected from X-ray point 
sources in RASS, and since clusters of galaxies (which 



represent the most massive DMHs) are also X-ray 
sources, an observational bias that might affect our re- 
sults is that AGNs in the clusters are selected against 
in the RASS-AGN sample due to confusion, especially 
considering that the mean point spread function (PSF) 
has a FWHMw 2 arcminutes (La Barbera et al. 2009), 
which is not negligible when compared to the extent of 
the cluster X-ray emission. We have estimated the sig- 
nificance of this effect as follows. The detection limit of 
the cluster diffuse X-ray emission is higher than that of 
point sources. A typical flux limit for clusters in RASS is 
Sx ~ 2 x 10~ 12 erg s" 1 cur 2 , (0.1-2.4 keV) (Bohringer et 
al. 2001; Popesso et al. 2004), corresponding to log Lx > 
44.3 [/if 2 erg s _1 ] at the lower bound of our redshift range 
of z w 0.16. This scales to a cluster virial mass of 
log M h > 14.5ft _1 M Q (Reiprich & Bohringer 2002). As 
we found at the end of Section 4.1, our CCF is only sensi- 
tive to the HOD behavior at log M h < 14.0/i _1 M Q , and 
thus this selection effect should play negligible role in our 
results. Additionally, a confusion due to two or more X- 
ray AGN point sources being artifically blended into a 
single source by the RASS point spread function (PSF) 
is not important either, as (N\) is much less than unity 
over all masses (see Fig. 5). Thus we can state with a 
good degree of confidence that our limits on a (Model A) 
or a s (Model B/C) are not an artifact of X-ray source 
confusion. 

Our results of a s < 1, found in our extensive explo- 
rations in the parameter space of the models that we con- 
sider, may be interpreted in view of the long-suggested 
deficiency of emission line diagnostic-selected AGNs in 
rich clusters (e.g. Gisler 1978; Dressier et al. 1985). While 
our CCF analysis is not sensitive to rich clusters (Sect. 
4.1), a more recent work by Popesso & Biviano (2006) 
verified that this trend extends to the group scale by 
finding an anti-correlation between AGN fraction and 
the galaxy velocity dispersion in nearby (z < 0.08) 
groups/clusters of galaxies down to the velocity disper- 
sion of <7 V ~ 200 [kms -1 ]. Arnold et al. (2009) also found 
that the X-ray selected AGN (log L x > 41/i 70 1 crg s _1 , 
representing a much lower luminosity population than 
our AGNs) fraction is larger in groups than in clusters 
by a factor of two at low redshifts (0.02 < z < 0.06). 
Various physical mechanisms have been suggested to ex- 
plain the deficiency of star-formation galaxies in clusters 
of galaxies, and similar processes might be responsible 
for the lack of AGNs as well. These include ram-pressure 
stripping (Gunn & Gott 1972) or evaporation (Cowie & 
Songaila 1977) of a galaxy's interstellar medium, which 
is a required ingredient of the AGN activity in galaxies, 
by the hot intracluster medium. Also, if AGN activity is 
triggered by mergers between two gas-rich galaxies, an- 
other possible mechanism for this trend is the decreased 
cross-section for galaxy mergers in the galaxy-galaxy 
close encounters at high relative velocities (Makino & 
Hut 1997). In high mass DMHs, representing richer 
groups/clusters, the velocity dispersion is higher and 
thus one may expect that there are fewer mergers than 
in groups of galaxies (Popesso & Biviano 2006), although 
the density of galaxies increase at the center. Other 
mechanisms such as the accretion of small gas-rich dwarf 
galaxies (minor mergers) may play a role in fueling the 
AGN host, but it is not clear how this would affect the 
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slope a s . 

In order to fully understand the growth history of 
SMBHs as well as the physical processes responsible for 
the AGN activity in groups and clusters, we need to 
explore the statistics at different redshifts, luminosities, 
and AGN types. The results mentioned above that inves- 
tigate AGNs in groups/clusters are concerned with AGN 
at lower redshifts and with lower luminosities than the 
AGNs in our CCF analysis. In the AEGIS field, Geor- 
gakakis et al. (2008) found that AGNs at z « 1 are more 
frequently found in the group environment than the over- 
all galaxy population. However, deep field surveys such 
as AEGIS do not contain a large enough volume to ex- 
plore the richness dependence of the AGN fraction in 
groups/clusters. Martini et al. (2009) investigate Chan- 
dra images of a sample of 33 clusters at 0.05 < z < 1.3, 
finding that the luminous AGN (log Lx > 43/i7 1 erg s _1 ) 
fraction increases with redshift. They explore the red- 
shift dependence of the AGN fraction, but not the rich- 
ness dependence. At z < 0.4, where the luminosity range 
of their AGNs is closer to ours, they have only two AGNs 
in 17 clusters in the redshift range that is similar to 
our sample. A much larger sample of groups/clusters 
is needed to explore the richness dependence of the AGN 
fraction in a sufficiently narrow redshift range to avoid 
any degeneracy with the redshift dependence. 

In general, directly investigating AGN prevalence in 
individual groups and clusters is observationally and an- 
alytically intensive. The HOD analysis of galaxy- AGN 
CCF provides a strong statistical tool with which to 
probe the AGN population in groups and clusters, which 
can be made without identifying individual groups and 
clusters. Thus it is important to extend the HOD anal- 
ysis to AGN clustering measurements at higher redshift, 
especially in view that the fraction of star-forming galax- 
ies in clusters increase with redshift (Butcher & Oem- 
ler 1984). It is also important to have larger numbers 
of AGNs for the CCF measurements, which will en- 
able measurements of small-scale clustering and produce 
stronger constraints on a s , as well as to break the model 
degeneracy. Using non-LRG galaxies for the CCF to ob- 
tain the 1-halo term constraints in the lower Mh regime 
would also improve the analysis. 

Linked to our constraints on a in Model A, we also 
have weak constraints on the range of M CI . The range 
of M cr is similar to that of M m i n in Model B. Using the 
scaling relations and the luminosity range of our RASS- 
AGN samples, we can make an interesting comparison. 
Figure 1 of paper I shows that the minimum 0.1-2.4 keV 
luminosities of the total and high Lx RASS-AGN sam- 
ples are log Lx [/if 2 ergs _1 ] « 43.7 & 44.3 respectively, 
while the low Lx sample has the same minimum lumi- 
nosity as the total RASS-AGN sample. Converting the 
0.1-2.4 keV luminosity to a 0.5-2 keV luminosity, assum- 
ing a power-law spectrum with photon index T = 2.5 and 
applying the luminosity-dependent bolometric correction 
from Eq. (2) of Hopkins et al. (2007), the corresponding 
bolometric luminosities are log Lboi [^7o 2er § s ~ 1 ] ~ 44.5 
and 45.4 for the total and high L x RASS-AGN samples. 
The mass of the SMBH (M.) corresponding to the min- 
imum luminosity is then M. w 2 x 10 6 /A and 2 x 10 7 /A 
[M ], where A = iboi/^Edd and L^dd is the Eddington 
Luminosity. 



The relation between the black hole mass M. and the 
DMH of its host galaxy is explored by Ferrarese (2002) by 
comparing the rotation curves of spiral galaxies. While 
their M.-DMH mass relation is subject to uncertain- 
ties due to different assumptions of the circular veloc- 
ity and virial velocity by a factor of a few, using their 
Eq. (6) with Bullock et al. (2001) 's recipe, the above 
black hole mass corresponds to a galaxy halo mass of 
4 x 10 11 (1.5 x 10 12 ) -A- O - 6 [M ] (for h = 0.7). These 
masses are comparable or slightly lower than the lower 
bounds of M cr derived above for those AGNs emitting at 
the Eddington luminosity and are an order of magnitude 
lower than our best-fit values. In general, the concept 
of the DMH of a galaxy explored by Ferrarese (2002) is 
different from the DMH derived in our CCF/HOD anal- 
ysis, since the large-scale correlation function reflects the 
DMH mass as the largest virialized structure the object 
belongs to. If the largest virialized structure represents a 
group or cluster of galaxies, then there are multiple ob- 
jects in a DMH, and the DMH that is related to Ferrarese 
(2002) 's scaling relation represents a sub-halo in a larger 
mass halo. While our statistics do not allow for strong 
constraints, if M cr is close to our nominal value, then 
the lowest luminosity AGNs among our samples should 
be accreting at a sub-Eddington rate (A < 0.1) or reside 
in groups. 

6. CONCLUSION 

We have applied a halo occupation distribution (HOD) 
model to the two-point cross-correlation function (CCF) 
between broad-line RASS-AGN and luminous red galax- 
ies (LRGs) at 0.16 < z < 0.36 measured in paper 
I. We have developed a method of applying the HOD 
model directly to the CCF, where we have used the well- 
determined HOD of the LRG as a reference to constrain 
the HOD of the AGNs. Our findings are summarized as 
follows: 

1. From our HOD analysis of the AGN-LRG CCF, 
we have found constraints on the HOD of AGNs. 
We have modeled the AGN HOD as a truncated 
power-law model (Model A) with two parameters: 
the critical dark matter halo (DMH) mass (M cr ) 
below which the AGN HOD is zero (i.e, there are 
no AGN in halos below this mass threshold), and 
the slope a, where the HOD oc for M h > M cr . 
We have obtained joint constraints in this two- 
parameter space, M cr — a under the scenario that 
there is no central AGNs in high mass halos. 

2. The constraints in this M cr — a space lie mainly 
along lines of constant large-scale bias and constant 
mean DMH mass. These provide a more accu- 
rate determination of the mean mass of the DMHs 
in which AGNs reside than the analysis based on 
power-law fits provided in Paper I. The mean DMH 
mass we derive is (M h >[/i _1 M Q ] = 13.09 ± 0.08, 
13.26±g;g|, and 12.97±g;^ for the total, high L x 
and low Lx RASS-AGN samples respectively. The 
2-halo term of the CCF primarily constrains the 
mean DMH mass. 

3. In all three of our RASS-AGN samples, we have de- 
tected a significant contribution of the 1-halo term 
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in the AGN-LRG CCF, indicating that some AGNs 
are in the same DMHs occupied by LRGs. This 
gives additional constraints on the AGN HOD be- 
yond the mean DMH mass. 

4. The range of acceptable HODs of high Lx and low 
Lx RASS-AGN samples are significantly different 
in the M CI — a space, confirming the paper I's re- 
sults on the X-ray luminosity dependence of AGN 
clustering properties at the ss 2a level. 

5. The most important constraint we have obtained 
using Model A is an upper limit (corresponding to 
Ax 2 = 2.3), of a < 0.4 for the total AGN sam- 
ple, rejecting a = 1. In the extreme case where 
all AGNs are satellites, the results show a sharp 
contrast with previously derived HODs of galax- 
ies, which show a « 1 — 1.2 (i.e., the number of 
satellite galaxies approximately oc Mh), across a 
range of galaxy luminosity and redshift. Taken 
at face value, this implies that the AGN fraction 
among satellite galaxies decreases with increasing 
Mh. This is consistent with previous observations 
that the AGN fraction is smaller in clusters than 
in groups in the nearby Universe. Possible expla- 
nations include the ram-pressure stripping and/or 
evaporation of galactic gas by the hot intracluster 
medium and a decreased cross-section for galaxy 
merging in the high velocity dispersion environ- 
ment of richer groups/clusters. 

6. We also investigated a model (Model B) which is 
composed of central and satellite AGNs, where the 
central HOD is constant (= /a) and the satel- 
lite HOD has a power-law form = / A (M/Mi) Qa , 
both at masses above M m - ln . This model is based 
on the results of HOD analyses of galaxies. For 
Mi/M m ; n = 23, which is appropriate for galaxies, 
we obtain a s < 0.95, finding a marginal prefer- 
ence for a picture in which the AGN fraction among 
satellite galaxies decrease with DMH mass. With 
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